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^ ; Abstract 

Q ■ In this work we study the hadron-quark phase transition matching relativistic hadro dynamical 

Ph . mean-field models (in the hadronic phase) with the more updated versions of the Polyakov-Nambu- 

Jona-Lasinio models (on the quark side). Systematic comparisons are performed showing that the 
■ predicted hadronic phases of the matching named as RMF-PNJL, are larger than the confined phase 

^SJ ■ obtained exclusively by the Polyakov quark models. This important result is due to the effect of 

the nuclear force that causes more resistance of hadronic matter to isothermal compressions. For 
' sake of comparison, we also obtain the matchings of the hadronic models with the MIT bag model, 

named as RMF-MIT, showing that it presents always larger hadron regions, while shows smaller 
^ i mixed phases than that obtained from the RMF-PNJL ones. Thus, studies of the confinement 

I transition in nuclear matter, done only with quark models, still need nuclear degrees of freedom to 

be more reliable in the whole T x fx phase diagram. 

PACS numbers: 12.38.Mh,25.75.Nq 



1 



I. INTRODUCTION 



The study of the properties of the strongly interacting matter is experimentally supported 
by the heavy-ion collisions at ultrarelativistic energies, accomplished in the most sophisti- 
cated accelerators such as the relativistic heavy-ion collider (RHIC) and the large hadron 
collider (LHC). The measurements coming from these experiments are indirectly used to 
furnish the basic information on the phase diagram concerning the state of matter, being 
the different regions related to the distinct hadronic (confined quarks), and the quark-gluon 
plasma (free quarks) phases. In order to cover the temperature-density/chemical potential 
plane of the hadron-quark phase diagram, the new experiments should be able to compress 
the baryonic matter even to higher densities compared to the well known nuclear saturation 
density, po ~ 0.16 fm~^. The new facility for antiproton and ion research (FAIR) [ij at 
Darmstadt, and the nuclotron-based ion collider facility (NICA) at the Joint Institute 
for Nuclear Research (JINR) in Dubna, will make possible such extreme conditions through 

I I 

collisions foreseen to reach an energy range of 8 to 45 A GeV [1]. Therefore, it will be possi- 
ble to test the predictions from the theoretical calculations about the order of hadron-quark 
phase transition (crossover, first order or both), and also the regime of high density and 
moderate temperatures. 

On the theoretical side, the strongly interacting matter is described by Quantum Chro- 
modynamics (QCD). In its nonperturbative regime of large distances, or equivalently low 
energies, the^ost important method to study the QCD structure is the numerical lattice 



calculations [3], with the Monte Carlo simulations being its mainly representative. Such 
techniques provide results as for the pure gluon sector, i.e., in the limit of infinitely heavy 
quarks, as for systems including dynamical quarks. The difficulties in the latter are in the 
fact that at finite quark chemical potential (fig) regime, the numerical calculations face to 
the fermion sign problem jsl. However, alternatives to solve this question are addressed in 
many methods On the other hand, one can deal with the high energy QCD regime 



(asymptotically free region) [Uj in a complete theoretical way through the perturbative 
treatment. 

To investigate the moderate and hi gh d ensity regime of the transition, it is needed effective 



models such as the MIT bag model 



12| and the Nambu-Jona-Lasinio (NJL) 13Nl5| one. 



that curiously were firstly proposed to describe the hadronic mass spectrum. The former 



2 



treats gluons and massless quarks as free particles in which the confinement phenomenon is 
incorporated by including a bag constant in an ad hoc fashion. The latter presents further 
similarities with the full QCD theory but do not take into account the confinement, since 
the quarks interact each other via pointlike interactions with no mediator gluons. One of 
the mainly aims of these kind of models is to provide the QCD phase diagram in the T x n 
plane in order to distinguish hadronic phase where the chiral symmetry is broken, and the 
quark one. The quark stars description with these effective models is in the same important 
way liliel. 

The phenomenology of the 
tioned effective QCD models 



ladron-quark phase transition is often studied with the men- 



liM aud other. Recently, a study about the QCD phase 

diagram was accomplished in Ref. [19| where the linear a model was used to construct the 
phase transition curve, in the two-fiavor quark system, as a function of the vacuum pion 
mass. It was shown that the point in the T x ^ plane where is changed the phase transition 
order, the critical end point, is sensitive to this parameter. In such cases, where there is 
only one model (the quark one), the identified phases in the transition curve are associated 
to those with broken and restored chiral symmetry. 

To distinguish hadron and quark phases we have the following well accepted phenomenol- 
ogy. The strongly interacting matter at very high temperatures and baryon densities should 
lave quarks and gluons as the degrees of freedom. Motivated QCD models, cited above 
12h15| take these degrees of freedom into account. On the other side, at low temperatures 
(T < 20 MeV) and moderate densities of the order of po, the hadronic phase can not be 
described by quarks and gluons. Instead, nucleons and mesons are the relevant degrees 
of freedom. Here, different approaches as Brueckner-Hartree-Fock by using realistic two- 
nucleon interactions 2^, or relativistic mean-field (RMF) models may describe quite well 
the nuclear matter and finite nuclei properties. Between these two phases a mixed region 
takes place. 

In this work, we focus in two structurally distinct models to treat the hadron-quark phase 
transition, namely, the Polyakov Nambu- Jona-Lasinio (PNJL) models in the characterization 
of the quark phase, and the relativistic mean-field nonlinear Walecka ones in the description 
of the phase where the chiral symmetry is broken, i.e., the hadronic phase. The former 
is a generalization of the NJL model in the sense that the confinement phenomenon is 
incorporated in its structure. For sake of comparison, we also perform a hadron-quark 
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phase transition using the MIT Bag model in the quark side, in order to see exphcitly the 
effect of the dynamical confinement contained in the PNJL model compared with the MIT 
Bag one. 

The motivation for such an approach matching two models with different degrees of 
freedom is the following. MIT bag and PNJL models do not describe any of the very well 
known properties of the ordinary nuclear matter while attempting only to achieve, through 
phase diagrams, the boundaries where this matter starts to be confined. In this aspect only 
with such models a more complete description is poor. On the other side, the RMF models 
constructed to well describe the bulk nuclear matter and finite nuclei properties, are often 
applied to investigate high density regimes as it is common for neutron stars. Whether in 
such regions the RMF models are still valid is questionable. The phase transition between 
the quark and the hadron phase helps understanding when one needs to start describing 
nuclear matter with other degrees of freedom than baryons and mesons employed by RMF 
models. Therefore, our study brings more information for the mixed phase between these two 
still disconnected quark and hadronic phases. The use of different PNJL and RMF models 
is needed if one wants to have a more precise conclusion about how the different families of 
hadronic and quark models predict mixed phases. If they are similar or not, or still whether 
the description of the confined/deconfined phase transition only via the PNJL model itself 
is totally trustable when compared with the approach matching the two models. As we 
will show, this phase transition obtained by only the PNJL model, is expected to occur in a 
much smaller chemical potential for the same temperature differing strongly of the transition 
obtained by RMF-PNJL models. Future experiments, discussed already before, will show 
what description is correct. 

After the first version of this manuscript was submitted, we became aware of a very 
similar study in which the hadron-quark phase transition was investigated only with one 
parametrization for the PNJL and RMF model 2l|]. Our systematic study here goes beyond 
that work, in the following aspects. Our calculations will take into account here different 
parametrizations for both, the PNJL and RMF models. In the latter, we will analyze a class 
of hard and soft equations of state, in a very large range of the nuclear matter incompress- 
ibility, 172 < K < 554 MeV, and show that, actually, the RMF models affect quite small the 
entire hadron-quark phase transition. In the former, the most up to date Polyakov potential 
parametrizations will be used and compared each other. Still, a comparative study involving 
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the hadron-quark phase transition with the MIT bag model will be performed. 

This paper is organized as follows. In Sec. II, we present the quark models used in this 
work, starting by the MIT Bag model and later with the PNJL one, where a comparison 
between the current parametrizations, regarding their thermodynamics, is also done. In Sec. 
Ill, the mainly features of the RMF models are presented, and in Sec. IV, the phase diagrams 
concerning the hadron-quark phase transitions are shown. Our summary and conclusions 
are in the last section. 



II. QUARK MODELS 
A. MIT Bag model 

Possibly, the simplest model to describe an approximate physics describing the matter 
where quark and gluons are the proper degrees of freedom of the system is the MIT bag model 



12| . Such a scenario is supposed to exist at very high density or temperature regimes. Back 
in the big bang theory, at the very beginning, the universe was very hot and dense before 
hadronization. Nowadays, very high energy heavy-ion collision experiments in LHC try to 
recreate the signatures to confirm this hypothesis. In the case it is confirmed, hot hadronic 
matter undergoes a phase transition to the quark-gluon plasma (QGP) at uncertain but very 
high and very low temperature and density, respectively. The phase diagram connecting the 
QGP phase to a pure hadronic phase from high to zero temperature is still a big challenge 
and has a model dependence exhibited in previous studies when some relativistic hadronic 



models and the MIT bag ones were investigated [2^, |23|]. The MIT bag model itself does not 
describe the confined hadronic phase but rather the QGP. Therefore, the hadronic phase 
has to be represented by any hadronic model in which the degree of freedom are baryons 
and mesons as we will present later on. Once each phase is modeled, the phase diagram is 
obtained thermodynamically by the Gibbs criteria, matching the chemical potential and the 
pressure at a given temperature (Tc), across the phase boundaries when a phase transition 
takes place. Explicitly, these criteria are 

fiH{Tc,PH) = f^QGp{Tc, Pqgp) (1) 
P^,(T„p^) = PQGp{n,p'QGp)- (2) 



5 



As previously presented 22|], the set of critical p^'s establishes curves p'^{T) in the {p,T) 
plane. Below the {p'j^,Tc) curve the system may be interpreted as nuclear matter described 
by the hadronic models (hadronic sector). The region in between (p^,Tc) and (pQCpyTc) 
corresponds to a mixed (H-QGP) coexistence phase. Above the {Pqgp^ ^c) curve, the system 
is in a pure QGP phase. 

For the quark-gluonplasma phase in the MIT bag model, the pressure and baryon number 
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density are given by 
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where B is the bag constant, Nf is the number of flavors and ag is the QCD running 
coupling constant, depending on the quark-gluon plasma temperature Tq and the quark 
chemical potential fiq through the first order perturbative expression 
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ln[(0.8p^ + 15.622r2) /A^ 



(5) 



In this paper, the MIT bag model will be used only as a comparison with the PNJL model, 
regarding the predictions for the boundaries of the hadronic and QGP phases for each specific 
RMF model describing the hadronic phase. 



B. PNJL model 

1. Confinement 

The confinement can be measured by the Polyakov loop, defined as a Wilson line in the 
Euclidean space-time by L = 'Pexp(i dTA4) where A4 is the gauge field. By this definition, 
it can be viewed as a closed path connecting the same point in the space, at the different 
times and /3 = 1/T. The traced Polyakov loop, $ = Ttc{L)/Nc with Nc being the quark 
color number, plays an important role concerning the deconfinement since it is the order 
parameter of this transition. Actually, the deconfinement is a consequence of the global 
center symmetry breaking. In pure gauge systems at finite temperature, where the periodic 
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boundary conditions are satisfied, this symmetry is realized by the gauge transformations 
due to the way how they affect the boson fields. The consequence in the traced Polyakov loop 
is that $ — )■ e $j with k = 1,2, 3... what means that the condition $ = has to be fulfilled 
in order to keep the center symmetry preserved. Since $ is related with free energy of a color 
source, through $ = e^^''^'^'^'^ , if the center symmetry is realized, then Fg — )■ oo, meaning 
that the quark is confined. In dynamical quark systems where the fermionic field explicitly 
breaks the center symmetry, generating automatically $ 7^ 0, the traced Polyakov loop is 
considered as an approximated deconfinement order parameter. The situation is analog to 
the cases in which the system Lagrangian density presents a mass term that explicitly breaks 
the chiral symmetry. The condensate (%l)qip^ is also an approximated order parameter, in 
this case for the restored chiral symmetry phase. The confinement effect, associated with 
$, was implemented originally in the NJL model by Fukushima and deeply studied in 



Refs. 



27 



33| . The PNJL model is able to describe as the broken- unbroken chiral symmetry 



as the confinement-deconfinement phase transitions even being this situations, respectively, 
in the opposite regimes of vanishing and infinite quark masses 26 1. 



2. Thermodynamical calculations at T ^ 

The connection between the fermion and the gauge field in the PNJL model is done by 
making the substitution + iA^ in the Lagrangian density. The after bosonization of 

the system, and the mean-field approximation lead to the following grand canonical potential 
per volume ^|: 

2 27T^ Jq 

r In [1 + 3<l>e-(^'-^')/^ + 3$*e-2(^«-'^«)/^ + e-3(i?,-/.,)/Tj ^2^^ 
Jo 

POO 

/ In [1 + 3<l>*e-(^^+'^^)/^ + 3<l>e-2(^^+'^^)/^ + e-=^(^^+'^^)/^] edk, (6) 
Jo 



2-k'^N, 
2n^N, 



being Eg = {k"^ + Mq^Y^"^ , psq = (^tpqipq) = {uu) + i^dd') = 2 {uu) in the isospin symmetric 
system, W($,$*,T) the effective Polyakov loop potential in terms of $ and its conjugate 
$*, that we will discuss later, and '~fq = Ng x Nf x = 12 the degeneracy factor due to 
the spin, flavor, and color numbers (A^^ = Nj = 2 and = 3). The constituent quark mass 
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should obey the autoconsistent equation 



M, = Mo- Gp,, (7) 
where the quark condensate p^g, determined by the condition {dVLp^^i^/ dpsq) = 0, is given by 

poo J. J /»A 

/ ^edk [F{k, T, $, $*) + F{k, T, $, $*)] - / -^edk (8) 



P. 



7g 



27r2 
with 



Fik T $ $*) = — — flOl 



being the generahzed Fermi-Dirac distributions. As pointed out in Refs. [34, 135|], an impor- 
tant consequence of the structure of the couphng between $ and the quark sector, is the 
possibihty to deal with the PNJL model in the same theoretical way as in the NJL one, 
regarding the statistical treatment. The modification in the equations, as, for example, in 
Eq. (IH]), is in the use of these new distribution functions for quarks and antiquarks. 

Through fipNjL ? all the thermodynamical quantities can be obtained, namely, the pres- 
sure, given hj P = —VL that leads to 

PpN.. = -w($,$*,r)-%! + -^ l\e + MlYi^kHk 

■y r°° k"^ 
an 

dfi ' 



the quark density, p = 



and the energy density, S = — T^^^S^ + pp 



2 /-A 



/■oo 



7. 



Notice here that we have subtracted the vacuum value of fipNjL, resulting in the addition of 
the VL^^^ term in the expressions (fTTl) and ([T3l) . 
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The entropy density is obtained through 5 = — |p, or by the thermodynamical relation- 
ship S = {P + £ — ^p) /T. Therefore we have 
f)jj ^ k'^ 



9r 67t-T j„ (B + w^yft 

POO 

/ + m2)1/2 k'^dk [F{k, T, $, $*) + F{k, T, $, 
Jo 

POO 

/ k^dk[F{k,T,^,^*) - F{k,T,^,^*)] . (14) 
Jo 



27r2T 



27T^T 

It should be mention at this point that the prescription of generalization of the NJL 
model to the PNJL one concerning the Fermi-Dirac distributions functions of parti- 
cles, f{k,T), and antiparticles, f{k,T), to those shown in Eqs. (jH]) and ffTU]) . is vahd 
for the entropy density functional given in Eq. (1141) . and not for its usual expression 
S ^Jd^flnf + (1 - /)ln(l - /)] + /In/ + (1 - /)ln(l - /)]. 

To deal with the entire PNJL model, characterized by the above equations of state (EOS), 
it is needed specify the Polyakov loop potential, W($,$*,T). Different versions were pro- 



posed in the literature, and following the language of Ref. 



RTW05 [2711, and RRW06 



2^35 



we refer two of them b\ 



35| . We call the other two used in our work by FUKU08 



and DSIO '37|. Their functional forms are given, respectively, by, 

2 -<I)<1>* + 64(r)ln [l -6$$* +4(<l)3 + $*^) -3($<1>*)2] , (16) 

= -546-"/^$$* + In [1 - 6$$* + 4($3 + _ 3($$*)2] , (17) 

Wdsio = (aoT^ + a.fil + a^T^^D^^ + asT^^ln [l - + 8$^ _ 3$^] , (18) 



y4 2 6 ^ ^ 4 



where 



h{T) = ao + a,\^-^j+a,l^-^j + «3 (^^ j , and b,{T) = b, 1^-^ j . (19) 

In a general way, the Polyakov potentials are constructed in order to reproduce the 
well established data from lattice calculations of the pure gauge sector (where $ = $*), 
concerning the temperature dependence of the traced Polyakov loop and its first order 
phase transition, characterized by the jump of $ from the vanishing to a finite value at 



To = 270 MeV (see the dotted curve of Fig. 2 in Ref. [28|, for instance). Actually, in 



this work we reduced this value to Tq = 190 MeV, following the rescaling adopted in Ref. 



27l |. in order to better reproduce the lattice QCD results that we will show in Fig. [21 and 
the transition temperature at vanishing chemical potential, Tc(/i = 0), also predicted in 
the lattice, and that we will discuss in Sec. IIVI By taking this rescaling into account, we 
match more realistic PNJL models with the RMF hadronic ones, that already describe quite 
realistically nuclear matter and finite nuclei. 

The FUKU08 potential presents two free parameters, a = 664 MeV and b = 0.03A^, and 
was derived from strong coupling lattice expansion. Also to correctly reproduce the lattice 
results for Tc(/i = 0) we changed here the b parameter value to 6 = O.OOTA'^. 

The RTW05 parametrization was based on a Ginzburg-Landau ansatz, presenting the 
polynomial form in terms of the order parameters, $, and $*. The improved version RRW06 
uses, instead the polynomial terms of third and fourth order, the logarithm of the Jacobi 
determinant. Differently from the FUKU08 potential, the dimensionless parameters a^, and 
bi of the RTW05, and RRW06 ones (see their values in Table [T|) are found with the additional 
fit to the lattice data for the energy density, entropy density, and pressure of the gauge sector 
as a function of the temperature, including the proper Stefan-Boltzmann (SB) limit at high 



temperature regime. Fig. 1 of Refs. [27| and [28|], for example, shows the good agreement 
among these quantities and the lattice results. 



TABLE I: Dimensionless parameters of the potentials given in Eqs. (jl5p -(|16p and (jlSp . 



Potentials 


ao 


ai 


a2 


«3 


63 


h 


RTW05 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 


RRW06 


3.51 


-2.47 


15.22 






-1.75 


DSIO 


-1.85 


-1.44 X 10-^ 


-0.08 


-0.40 







The DSIO potential, proposed in Ref. [37|, presents baryon chemical potential dependence 
and is used in a hybrid SU{3) chiral model that has both hadrons and quarks as degrees 
of freedom. This parametrization is able to reproduce the gauge lattice results, just like in 
the potentials shown before, and also allows that the hadron-quark phase transition occurs, 
at zero temperature, at density of 4 times the saturation density, and that the critical 
end point (CEP) is situated at jjB = 354 MeV, and T = 167 MeV By using this 

particular parametrization, we remark that the quark density, Eq. ( IT2|) . and the energy 
density, Eq. ( IT3l) . have to be modified, respectively, to Pq ^ Pq — (4ai/Xq + 2a2T^/Xq)$^ and 
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^PNJL ^PNjL - (4ai/i^ + 2a2TVg)$^ • 

To define completely the PNJL model it is needed to determine the coupling constant 
G, the cutoff A, and the current quark mass Mq = = M^. This is done by imposing 
the reproduction of the certain vacuum values, namely, the pion decay constant fixed in /^r, 
the quark condensate (uu), and the pion mass m^r. These values, together with the zero 
temperature expressions of the PNJL model, are shown in Appendix A. 



3. Comparing the Different Polyakov Potentials 

To evaluate the PNJL model one needs to solve simultaneously Eq. ([7]) and the mini- 
mization conditions for the thermodynamical potential relatively to $ and $*. Along all our 
study we will follow the lowest order approximation described in Refs. 30|, l32| that auto- 
matically leads to $ = $*. This approach reduces the set of coupled equations to Eq. (JTj), 
and 



dU{^,T) 3T7, 



9$ 27r2iV, 



/ edk[g{k, T, $, /ij + g{k, T, $, -/i,)] = (20) 

^0 



with 



1 + (^-{Eq-^lq)/T 

coming from the condition (Sfip^jL 

/d<l>) = 0. 

To see how sensitive is the PNJL model regarding the different potentials, we present the 
behavior of some important quantities. First we show in Fig. [T]how the order parameters 
$ and psq vary as functions of the temperature for some fixed values of the quark chemical 
potential. 
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FIG. 1: ^ (la-lc) and the ratio psq{T) / psq{vac) (Id-lf), of the different PNJL parametrizations, 
as a function of the temperature for some fixed values of pq. 



If we use only the PNJL model to construct the T x fig diagram, the maxima o: 



and (dpsq/dT) at a given pq are used to define the transition temperature [27|, |30|, l32|, in 
this case a crossover transition. Indeed, the strength of this mixing is utilized to find the 
value of b in the FUKU08 potential 36|. Other similar criterium to find the transition point 
is to localize the maxima of the chiral, and Polyakov susceptibilities 26|, l32]. Here we will 
also use the Gibbs criterium to find such transition points since we are dealing with two 
different models to treat the different phases. 

In Fig. [2] we show, at vanishing chemical potential, the behavior of some important 
thermodynamical quantities. The pressure, energy density and entropy density as a function 
of the temperature, and scaled by the respective Stefan-Boltzmann values Psb = ^^^^f^^; 
^SB = 3PsB, and Ssb = 4:Psb/T, are displayed in Figs. 2a-2c. In units of T"^, we furnish 
the called interaction measure given hj S ~ 3P. This quantity show us how the model 
deviates from the noninteracting massless quark system, since in such regime its value is 
vanishing. Finally, as a function of the energy density, we see the evolution of the ratio P/S 
that is closely related with the sound velocity through = ^ = ^^^aJ^^ + ^- In the SB 



d£ 



limit, one has c^g^ 



Psb/^sb — 1/3. 
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FIG. 2: (Color online) (a-c): Pressure, energy density, and entropy density, in units of its respective 
SB values, versus T. (d): Interaction measure, in units of T^, versus T. (e): Ratio P/£ versus 



£^/^. All calculations at fiq = 0. The three flavor lattice results were extracted from Refs. 



38 



As a reference the recent QCD lattice results with temporal extent A',- = 4, 6 [38j, and 
AT, = 8 for the 2 + 1 flavor system, in calculations with the improved fermion actions 
asqtad, and p4, are also shown in Fig. O Rigorous comparisons among lattice data and 
quantities such as the difference ^ ^^f^ — and the moments of the pressure, can be 



found, for example, in Refs [271, l28|, |40| , and [28|, l30l-l32| , respectively. 



Notice here the good agreement of the PNJL parametrizations and the lattice data, even 
the latter being originated from the 3 flavor system. We highlight RTW05 and RRW06 as 
being the models that better agree with the data. These results can be explained by the 
rescaling procedure adopted in the Tq value (from 270 MeV to 190 MeV). 



As mentioned in Ref. [36[, the FUKU08 model was not constructed to reproduce the 
Stefan-Boltzmann limit in high temperatures. This explains the deviation of the lattice 
data from temperatures higher than 200 MeV even after our rescaling in the b parameter of 
the model, from b = 0.03A3 to 6 = 0.007A3. 
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III. RELATIVISTIC MEAN-FIELD MODELS 



Different from the microscopic approach in which the nucleon-nucleon potentials are de- 
veloped to reproduce the well established data of the few-nucleon physics, and where the 
basic informations about the many-body system are extracted, e.g., via Brueckner-Hartree- 
Fock methods, the Quantum Hadrodynamics (QHD), based on local Lagrangian densities, 
uses the nuclear matter bulk properties observables at zero temperature to adjust its free pa- 
rameters and thereafter construct all the thermodynamics of the particular hadronic frame- 
work. One of the most used models coming from QHD is the nonlinear Walecka (or Boguta- 

n 

Bodmer) [41| model, that is given by the following renormalizable Lagrangian density 

A B 

- g^aipip - g^il^i^uj^t^ - -a^ - -a^, (22) 

with F^i, = duOJ^ — dfj^ojy, and the nucleon degree of freedom represented by the spinor ip. 
The responsible mesons for the attractive, and repulsive parts of the nuclear interaction are 
described in this formulation by the scalar (a), and vector (w^) neutral fields, respectively. 
The saturation in nuclear matter is understood in the QHD by the almost vanishing value 
oi = S + V aX the saturation density, with 5, and V being the Lorentz meson potentials.. 
M, and m^^ are, respectively, the masses of the nucleon, scalar, and vector mesons. 



In their model, Boguta and Bodmer [4l| considered the cubic and quartic self-interactions 



in the scalar field a, in order to improve the original Walecka model [42|]. In this version, 
the models are able to control through the fitting of their coupling constants, the values of 
the saturation density, binding energy {Bq) as well as the incompressibility {K), and the 
effective nucleon mass (M*). 

Through the Dirac equation of the model, one identifies its nucleon effective mass as 

M* = M + g^a = M - Gl\ps + a{^Mf + h{AM)\ (23) 

with Gl = gl/ml, a = A/gl, b = B/g^, and AM = M* — M. Let us remark here that this 
definition of effective mass, also called Dirac mass, is valid for the relativistic models. A 
deep discussion about other definitions and concepts of this physical quantity can be found 
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in Ref. 43|. The scalar and vector densities are written, at finite temperature, as 

7 M* 



Ps 



p 



(^7V> 



27r2 

0_,A 7 



-^k'[f{k,T) + fik,T)]dk and 



e[fik,T)-fik,T)]dk, 



(24) 
(25) 



with 7 = 4 for symmetric matter. The usual Fermi-Dirac distributions to particles 



and antiparticles are defined, respectively, by f{k,T) 



,iy/k^+M*'^-u)/T 



+ 1 



and 



, and the relation between the effective chemical poten- 
tial, u, and the baryon chemical potential is given by = z/ + C^p. 

From the energy-momentum tensor obtained through Eq. (122]) . one can obtain the pres- 
sure of the system, that reads 

Glp^ (AM)2 a(AM)3 b{AMy 7 r k^ 



2GI 



6vr2 Jq (P + M*2)i/2 



p^^. = ^ — — ^^=^ — "'~'r' +^ I ———^T-^[f{k,T)+f{k,T)]dk, 

(26) 

already written in terms of the mean-field approximation, as well as Eqs fl23|) - fl25|) . 

For the construction of the hadron-qua rk p hase diagrams, we choose the following set of 



RMF parametrizations: Walecka (WAL 
NL3 



4^, NLB y 



NL2 j46|, NLSH NLBl j46|, 
. The WAL and NLZ2 models are 



47|, NLB2 |46|, NLC j45|, NLl [46| and NLZ2 |48 
chosen because their incompressibilities, K = 554.4 MeV and K = 172.0 MeV, respectively, 
indicate very hard and very soft EOS for infinite nuclear matter. The other models inter- 
mediate both. A full list of the nuclear matter saturation properties including po, Bq and 
m* = M*/M is found in Table HTl of Appendix B. 

Important features can be highlighted concerning these relativistic models. At moderate 
temperatures, T < 20 MeV, a liquid-gas phase transition is predicted In Figs. [3^-|3t 
we show the typical van der Waals behavior of the used parametrizations in the pressure 
versus density curves. 
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FIG. 3: (Color online) (a)-(c) Pressure as a function of p/po at T = 10 MeV, and (d)-(f) Effective 
nucleon mass as a function of the temperature at /u^ = 0. 



The Maxwell construction can be used to find the coexistence points, and consequently 
the coexistence curve. It was shown in Ref. jsol that this curve, scaled by the critical 
parameters, is not universal in the liquid region due importance of the interactions in such 
phase. Another interesting feature of the RMF models occurs at high temperature regime 
(T ~ 200 MeV), where other kind of phase transition takes place. This scenario was studied 
by Theis et al 511] in the context of the linear Walecka model at p = = 0. The authors 
have shown that the order of the phase transition depends on the values of the coupling 
constants used to fit po; and Bq. A signature of this transition is exhibited in the behavior 
of the nucleon effective mass as a function of the temperature, since an abrupt decreasing 
of M* characterizes a first order phase transition. In Figs. |3ji-|3]F we present the high 
temperature regime of RMF models used in our work concerning the behavior of M* at zero 
baryon chemical potential. A recent study regarding the high temperature regime of RMF 



models at different number of nucleons and antinucleons (p^ 7^ 0) was done in Ref. 



44 



51. 



IV. RESULTS AND DISCUSSIONS 



Now we present our mainly results regarding the hadron-quark phase diagrams, con- 
structed from the Gibbs criteria shown in Eqs ([I])- ([2]), defining first our notation for the 
transitions curves we will exhibit. We call the matching between the RMF models and the 
PNJL ones as the "RMF-PNJL" transition curves. When the transition is done with the 
MIT Bag model on the quark side, the modification is straightforward in the sense that the 
curve is denoted by "RMF-MIT" . Thus, the Gibbs criteria applied, for example, in the NL3 
model together with the RRW06 one, generate the curve called NL3-RRW06. If the PNJL 
model is replaced by the MIT Bag one in this case, we denote this transition by NL3-MIT. 
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A. RMF-PNJL transitions 



Regarding the condition of equal pressures, we stress that the hadronic pressure, Eq. (126|) . 
was used to match the quark ones, Eq. ([II]), and Eq. for the RMF-PNJL and RMF-MIT 
curves, respectively. 

Let us remark here that for the RMF-PNJL transitions, we adopt the Gibbs criteria with 
the additional constraint that the RMF confined phase is actually matching the deconfined 
phase of the PNJL model. We only consider the solutions satisfying this condition. An 
example of this restriction is shown in Figs. H^-Hb for the NL3-RRW06, and NL3-FUKU08 
matchings at = = 660 MeV. 
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FIG. 4: (a) Pressure as a function of the temperature for two of the PNJL models and the NL3 
one. (b): d^/dT and —d^/dT versus T, where ^ = psqiT) / psq{vac). The transition temperature 
(see text) is given by T. Both figures at fis = 660 MeV. 



Notice that the FUKU08 and NL3 pressures cross each other at the temperature given 
by T = 128.8 MeV, but since at this temperature the FUKU08 model is still in the confined 
phase (see the down panel of Fig. Hb), we did not consider this crossing for the NL3-FUKU08 
phase diagram. Here we are following the same procedure used in Ref. j27j concerning the 
definition of the transition temperature. In the case of the adopted rescaling for Tq (in 
RTW05, RRW06 and DSIO) and for b (in FUKU08), the almost perfect coincidence between 
the peaks of the temperatures derivatives of $ and psq is lost. To circumvent this problem, 
for each chemical potential, the transition temperature T, as the average of the different 
temperatures related to the peaks of d^/dT and dpsq/dT. For the FUKU08 potential at 
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fiB = 660 MeV, this average is T = 139.5 MeV. 

For the RRW06 and NL3 models, its pressures cross at T = 149.2 MeV, where the RRW06 
model already predicts the deconfinement, since its transition temperature is given by T = 
148.8 MeV (see the up panel of Fig.. |4)d). Therefore, the point in which T = 149.2 MeV 
and /i_B = 660 MeV contributes to the NL3-RRW06 phase diagram. 

In the next figure we present the hadron-quark phase diagrams for all the RMF-PNJL 
parametrizations used in this work. 
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FIG. 5: Hadron-quark phase diagrams in the T x hb plane for the RMF-PNJL matching (grey 
band), and for the PNJL models themselves (full -|- dotted lines). 



In Fig. [5l each individual panel depicts the diagrams where the matching between the 
models is done by keeping fixed the indicated quark model, and changing the hadronic one. 
In the same figure, for comparison, we also give the transition line of the phase diagram for 
each PNJL model itself constructed, as aforementioned, from the average of the temperatures 
related to the peaks of d^/dT and dpsq/dT. The fuH (dotted) lines stand for the crossover 
(first order) transitions. It is important to stress here that there are different criteria to 
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construct these PNJL transition curves, such as the choice between different peaks of d^/dT 



53| and dpsq/dT. This fact can favor the emergence of the quarkyonic phase, predicted in 



Ref. j54|. A detailed study about these different choices and its consequences in the hadron- 
quark phase diagram is under progress. 

An important feature regarding the PNJL curves are the values of the transition tem- 
peratures at vanishing chemical potential, Tc(/i = 0) = 179 MeV, 181.8 MeV, 177.3 MeV, 
and 174.5 MeV, respectively, for the RRW06, RTW05, FUKU08 and DSIO parametrizations. 
Notice that these values agree very well with the lattice QCD result for this quantity given by 
Tc(/x = 0) ~ 170 — 190 MeV jsS], or even with the more stringent value of Tc{p = 0) = 173±8 



MeV 



56| . This nice agreement is completely destroyed if the original values of Tq and h are 



used in the construction of the PNJL curves. In this case, the values of Tc(/i = 0) are higher 
than 200 MeV for all the Polyakov potentials. 

We also remark that the PNJL phase diagrams could also furnish other phases beside 
the confined/broken chiral symmetry, and deconfined/restored chiral symmetry ones, delim- 
ited by the full and dotted lines in Fig. [5l These specific phases are closely related to the 
instabilities of the ground state, present in any fermionic system at sufficiently low temper- 
atures, that are overcame by the formation of the fermion pairs (in BCS theory, these are 
the so-called Cooper pairs). In particular, the quark system described by the PNJL models 
should also be affected in this regime, giving rise to the emergence of the two-fiavor su- 
perconducting color (2SC), and the color flavor locked (CFL) phases. This phenomenology 
can be incorporated in the PNJL model, via inclusion of the diquark condensate term in its 
Lagrangian density. However, in this work we focus exclusively in the nonsuperconducting 
phases of the hadron-quark diagrams by treating only the simplest version of the PNJL 
model, in which the color condensates are not being taken into account. Very interesting 
studies concerning the consequences o 



the 2SC and CFL phases in the description of quark 
stars were performed for the NJL j57| model, in which it is shown the role played by the 
diquark coupling stren gth s. Recently, this effect of the color superconductivity was also 
analyzed for the PNJL [58] model, where a modified version of the DSIO Polyakov potential 
is used in the description of the quark core of massive hybrid stars. 

Regarding the RMF-PNJL transitions themselves, we stress the following interesting re- 
sults. First of all, notice that for all the RMF-PNJL matchings, one can delimit very narrow 
bands that encompass all the transition curves, being the WAL-PNJL and NLBl-PNJL 
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matchings their extremes. This is not an obvious result since we are deahng with a large set 
of RMF models, in the sense that they present a range of very soft to very hard equa- 
tions of state, i. e., from the NLZ2 model with K = 172 MeV, up to the WAL one 
in which K = 554 MeV. This result indicates that at least for RMF models with scalar 
self-interactions of third and fourth order, all of them will predict very similar behavior 
concerning the hadron-quark phase transition, in connection with PNJL models. 

Our systematic study still shows that the additional constraint exemplified in Figs. 
andHb, are useful to define the CEP's of the RMF-PNJL transitions, given approximately 
by (/i^^^P = 650 MeV, T^^^ = 149 MeV), (750 MeV, 142 MeV), (750 MeV, 126 MeV), and 
(750 MeV, 138 MeV), respectively, for the RMF-RRW06, RMF-RTW05, RMF-FUKU08, 
and RMF-DSIO transitions. In a different way, the CEP is also defined in Ref. |21|. Notice 
how this restriction avoids the hadron-quark phase transition curves go inside the confined 
phases predicted by the PNJL transitions. 

This situation is completely modified if we use the original values of Tq and b of the PNJL 
models. In this case there is no CEP for the RMF-PNJL transitions, in the sense that there 
is no transition inside the hadron phase of the PNJL models themselves. The only exception 
is for the RMF-FUKU08 transitions, that present the same qualitative behavior shown in 
Fig. [5] for the calculations done with the original value b = 0.03A^. 

We still highlight here a very important physical consequence of the construction of the 
hadron-quark phase transitions via connection of the RMF and PNJL models. The panels in 
Fig. [5] clearly show us that the hadronic degrees of freedom present in the RMF models, make 
the system much more resistant to the quark liberation, in function of the density /chemical 
potential increasing, when compared with a system only described by the PNJL model. The 
physical origin of this important difference is due to the very known repulsive nature of the 
nuclear force, represented in the RMF models by the vector interaction, which strength is 
controlled by the coupling constant g^}, see Eq. ( l22l) (in nuclear matter equations of state, 
this strength is actually controlled by the ratio = Quj/itii^). The repulsion makes the 
system support more strongly isothermal compressions. In other words, the RMF-PNJL 
transitions predict confined (hadron) phases larger than those obtained exclusively with the 
PNJL models. This region of highly compressed matter of the hadron-quark phase diagram, 
is expected to be reached in the new experiments, such as the planned to occur in the 
new facihties FAIR/GSI (l|, and NICA/JINR Q. Therefore, it will be explicit what is the 
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magnitude of the role played by the repulsive interaction part of the nuclear force, even 
guiding possible selections of better parametrizations. 

The difference between the RMF-PNJL and the PNJL transitions can be decreased if we 
use PNJL models that also contain repulsive interactions, i.e., that present vector fields in 
its structure. Notice that the PNJL models used here are based on a structure that present 
only attractive interactions. There is no explicit terms proportional to the quark density in 
Eq. ([6]), coming from vector- type fields. Moreover, as pointed out in Ref. 36|, there is no 
constraint at all for the choice of the strength of this kind of interaction in the PNJL model. 
A study about the determination of this magnitude, by making minimum the difference 
shown in Fig. El is underway. 

B. Comparison with RMF-MIT transitions 

We perform now, systematic comparisons between the RMF-PNJL transitions with 
RMF-MIT ones, in order to verify explicitly the role played by the dynamical confinement 
of the PNJL models in the hadron-quark phase diagrams. To construct such curves, we use 
the particular value of 5^/^ = 238 MeV, that furnish a critical temperature at vanishing 
chemical potential consistent with the lattice simulation results. 

Firstly, we show in Fig. Elthe behavior of RMF-MIT curves in the T x fi^ plane for all 
the hadronic parametrizations. 

Some interesting points about these particular transitions have to be mentioned here. 
Notice that differently from the RMF-PNJL diagrams, all the RMF-MIT curves start at the 
same temperature, Tc{fi = 0) = 168 MeV. There is no critical end points as shown in Fig. 
|5l Similar behavior between the curves can be seen at T = 0, since all used models are lied 
in a narrow band in this region. They lie in the range around 1190 MeV < fiB < 1310 MeV 
(1730 MeV < fiB < 1810 MeV) for the RMF-PNJL (RMF-MIT) curves. For the RMF-MIT 
transitions, this behavior is strongly changed if the RMF models with higher order terms in 
the vector field, or even mixing terms between a and Un are used on the hadronic side. This 



is the case, for example, for the models used in Ref. [22j. We stress that the behavior shown 
in Figs. [5] and [6] is characteristic of the transitions in which the used RMF parametrizations 
only contain cubic and quartic self-coupling terms in the scalar field a, see Eq. 



This almost model independent result for the hadron-quark transition at T = for the 
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FIG. 6: (Color online) Phase diagrams in the T x /x^ plane for the RMF-MIT matching. In the 
Eq. daD, we have used B^/^ = 238 MeV and = 0. 

RMF-PNJL/MIT matchings, see Figs. OandEl is not surprising from the point of view of the 
PNJL models, since the Polyakov potentials related to the RRW06, RTW05, and FUKU08 
models vanish in the zero temperature regime, see Eqs. (|T5 l) -( |T7l) . and the modified Fermi- 
Dirac distributions, Eqs. (|9]) and ( !T0|) . behave as the conventional ones in the T — )■ limit. 
In particular, the DSIO Polyakov potential vanish at T = also due to $ = 0, see Eq.. (fTSll . 
That is, all the PNJL models used here are converted in the same model at T — )■ 0, in this 
case, the NJL one on the quark side. This is also the case in the RMF-MIT transitions, i.e., 
one has only one parametrization of the quark model since we fixed B and in Eq. ([3]) 
to construct the diagrams shown in Fig. [6l We remark that by construction, at T = all 
the RMF models lead to different bulk nuclear matter properties. We stress here that this 
almost model independence in the fiB value at T = can be relevant to the study of hybrid 
stars, composed by a quark core and a hadron crust 59|, iGOj. 

In order to see the effect of the dynamical confinement present in the PNJL models 
in the RMF-PNJL transition curves, we explicitly compare the phase diagrams of both, 
RMF-PNJL/RMF-MIT phase diagrams in the same Fig. H 
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FIG. 7: RMF-PNJL/RMF-MIT phase diagrams in the Txfi^ plane for some hadronic parametriza- 
tions. Here we also use ag = 0. 

In this figure one can clearly see that the hadronic region predicted for the RMF-PNJL 
matchings is always smaller than that obtained from the RMF-MIT ones, i.e., the quark 
degrees of freedom emerge before in the RMF-PNJL curves. Although we have presented 
only the transition curves shown in Fig. [71 we streamline that all the other diagrams follow 
the same pattern, thus the respective RMF-PNJL/RMF-MIT curves can be considered as 
representatives of all the models treated in this work. An enlarged view of these features 
can be viewed in Fig. |8l 

In this figure, despite present the larger hadron regions, that are delimited by the left 
branches of the curves, the RMF-MIT diagrams show that the mixed phase containing both 
hadrons and quarks is actually smaller in such diagrams. This difference, specifically at 
T = is extremely important to the study of hybrid stars, since its quark core is directly 
affected by the maximum density of the mixed phase 59|, l60 |. 

As our final remark, we point out that in the construction of the RMF-PNJL diagrams, 
we proceed in a different way that used in Ref. 6l|. In this reference, the authors used a 
bag constant also in the PNJL pressure in order to ensure that the quark pressure is less 
than the hadronic one in the confined phase. Here we force Ppnjl(7' = = 0) = 0, by 
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FIG. 8: RMF-PNJL/RMF-MIT phase diagrams in the T x p plane. 



subtracting from the grand-potential its vacuum value. 



V. SUMMARY AND CONCLUSIONS 



In this work we have studied the hadron-quark phase diagrams using two different models 
in the description of the distinct phases. The quark degrees of freedom were described via the 
recently suggested PNJL model 26| which incorporates the confinement phenomenon in the 
previous NJL one. As one of our results, we compared the current PNJL parametrizations 



that differ each other by the Polyakov potential $*,T^ 



RTW05 |27|, the logarithmic RRW06 |28|-|35|, the FUKU08 



We used the polynomial form 
and also included the DSIO 



371] one that presents a chemical potential dependence. The Fig. |2]show the good agreement 



with the lattice data, specially for the RRW06 and RTW05 parametrizations even the PNJL 
models being treated in the two-flavor system. It was also shown that the parametrizations 
furnish similar results for the analyzed thermodynamical quantities. The good agreement 
with the lattice data remains valid even for the transition temperature at vanishing chemical 
potential, predicted to be given by Tc(/i = 0) = 173 ± 8 MeV js^. Our calculations give the 
values of T^{fi = 0) = 179 MeV, 181.8 MeV, 177.3 MeV, and 174.5 MeV, respectively, for 
the RRW06, RTW05, FUKU08 and DSIO PNJL models used in this work. We still remark 
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that these nice values are obtained with the rescahng of the parameters Tq (RRW06, RTW05 
and DSIO) and b (FUKU08). 

On the hadronic side, we have used the well known RMF nonlinear Walecka models in 
its version that contains cubic and quartic self-coupling in the scalar field. We chose to deal 
with a set of these models that encompasses several incompressibilities values, representing 
very hard and soft equations of state. 

Regarding the phase diagrams, we have constructed the RMF-PNJL curves using the 
Gibbs criteria and the additional constraint that the RMF pressure matches the PNJL 
one only in its deconfined phase. This condition is exemplified in Figs. S^-IDd, and 
played an important role in the final diagrams, since it determines the critical end 
points given approximately by (/ig^^ = 650 MeV, T^^^ = 149 MeV), (750 MeV, 142 MeV), 
(750 MeV, 126 MeV), and (750 MeV, 138 MeV), respectively, for the RMF-RRW06, 
RMF-RTW05, RMF-FUKU08, and RMF-DSIO transitions. Moreover, it is also important 
to stress that all these transitions furnish very similar results in a such way that one can 
define a very narrow band in the T x fi plane encompassing all the RMF-PNJL curves, being 
the WAL-PNJL and NLBl-PNJL the limiting curves of these bands, see Fig. |5l This is a 
surprising result since we are dealing with a very large class of RMF models. In principle, 
there is no reason to the hadron-quark phase diagram behaves in a very similar way with 
such variety of RMF parametrizations. 

Other important result concerning the RMF-PNJL phase diagrams is their different pre- 
dictions compared to those obtained exclusively with the PNJL quark models. We found 
that the hadron phase described by the RMF-PNJL transitions is meaningfully larger than 
that predicted by the PNJL ones, i.e., the RMF-PNJL hadron phase is more resistant to 
the isothermal compressions. Physically, such difference is due to the repulsive part of the 
nuclear force described in the RMF models by the vector field interaction. Therefore, one 
become clear that there are very different results in treating the hadron-quark phase transi- 
tion via two distinct models, taking into account the different degrees of freedoms (hadrons 
and quarks), or only via quark models, even being the latter very realistic ones as the PNJL 
models that nicely agree with lattice QCD data, and where the dynamical confinement is 
considered in the NJL model through the Polyakov loop. We also remark that the region 
where the different descriptions do not agree each other will can be accessed in the future 
experiments planned to occur in the new facilities FAIR/GSI [l|, and NICA/JINR j2|. 
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As a final result, we compared the RMF-PNJL curves with that constructed by the match- 
ing between the hadronic models and the MIT Bag one, that incorporates the confinement 
via inclusion of the bag constant B. The Figs. [7] and [H] show that the dynamical confine- 
ment predicted in the PNJL model force the hadronic phase of the RMF-PNJL diagrams be 
smaller than the RMF-MIT ones. Curiously, the opposite occurs for the mixed phase, where 
hadrons and quarks coexist, see Fig. |8l These comparisons were done by assuming the value 
of B^^^ = 238 MeV for the MIT bag model, that nicely predicts a value of 168 MeV for the 
transition temperature at vanishing chemical potential, see Fig. [61 
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Appendix A: Zero temperature expressions for the PNJL model: NJL sector 



In this appendix we show the zero temperature expressions of the PNJL model, that 
actually are exactly the NJL ones, and give the remaining parameters needed to define the 
PNJL models presented in our work. We refer here only to that parametrizations that do 
not present contributions containing the traced Polyakov loop in the T = regime, i.e., we 
only consider the cases in which $*, T = 0) = 0. 

In the T = regime, the energy density given in Eq. ( fT3l) . and the scalar density in 
Eq. (IE]), are replaced by 
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where kpq is the Fermi momentum of the quark. The quark density is pq = ^kp^, and the 
pressure reads 

Pp,,,,l(T = 0) = figpg - £p^,^{T = 0), (A5) 

with the quark chemical potential given by pg = {kpg + M^Y^"^. 
Thus, the respective vacuum expressions obtained at kpq = are 



= and (A6) 



A 



vac 



So, fixing the values = 140.51 MeV, = 94.04 MeV, and | {uu) l^/^ = 251.32 MeV, 
and using Eq. ( 1A7|) together with the Gell-Mann-Oakes-Renner relation mlf^ = —MoPsq, 
and 

, _ AT.Ar^M^^ edk 

one obtains A = 651 MeV, Mq = 5.5 MeV, and G = 10.08 GeV"^. The constituent vacuum 
quark mass obtained from these values is = 325.53 MeV. 

The assumption Sp^^ji^iT = pg = 0) = leads to the following final expression for the 
energy density. 



^ac, (A9) 



where f^^ac = ^^vac- Notice that the same condition also ensures that Ppnjl(T = pg = 0) = 0. 
For the parameters aforementioned one has that = 409.15 MeV. In this work we have 

considered this value of fi^ac in the Polyakov potentials, even for the DSIO parametrization. 



Appendix B: Parametrizations of the RMF models 

Some important saturation quantities of the RMF models used in our work are listed in 
the next table. 

The binding energy is calculated from the energy density at T = 0, 
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TABLE II: Saturation properties of the RMF models. 



Model 



pa (fm 3) 



Bo (MeV) 



K (MeV) 



Walecka (WAL) 

NLB 

NL2 

NLSH 

NLBl 

NL3 

NLB2 

NLC 

NLl 

NLZ2 



0.150 
0.148 
0.146 
0.146 
0.162 
0.148 
0.162 
0.148 
0.152 
0.151 



-15.75 
-15.75 
-17.03 
-16.35 
-15.74 
-16.30 
-15.73 
-15.75 
-16.42 
-16.07 



0.54 
0.61 
0.67 
0.60 
0.62 
0.60 
0.56 
0.63 
0.57 
0.58 



554.38 
420.00 
399.20 
355.36 
280.00 
271.76 
245.10 
225.00 
211.70 
172.00 



by Bq = Sjimf/p — M at p = pq. The incompressibihty, K = 9^, reads 



dp ' 



gClp + ^ 
E* 



^ + 2aAM + 36(AM)2 + 3 



Ps 

M* 



(B2) 



with E*p = {kl + M*2)i/2. 
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